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Recently we proposed a theory of point-contact spectroscopy and argued that the splitting of 
zero-bias conductance peak (ZBCP) in electron-doped cuprate superconductor point-contact spec- 
troscopy is due to the coexistence of antiferromagnetic (AF) and ei-wave superconducting orders 
[Phys. Rev. B 76, 220504(R) (2007)]. Here we extend the theory to study the tunneling in the 
ferromagnetic metal/electron-doped cuprate superconductor (FM/EDSC) junctions. In addition to 
the AF order, the effects of spin polarization, Fermi-wave vector mismatch (FWM) between the FM 
and EDSC regions, and effective barrier are investigated. It is shown that there exits midgap surface 
state (MSS) contribution to the conductance to which Andreev reflections are largely modified due 
to the interplay between the exchange field of ferromagnetic metal and the AF order in EDSC. Low- 
energy anomalous conductance enhancement can occur which could further test the existence of AF 
order in EDSC. Finally, we propose a more accurate formula in determining the spin polarization 
value in combination with the point-contact conductance data. 

PACS numbers: 74.20.-z, 74.25. Ha, 74.45.+C, 74.50.+r 



I. INTRODUCTION 



Using point contact technique to measure the spin 
polarization in ferromagnetic metal/conventional super- 
conductor (FM/CS) junctions was pioneeringly done by 
Soulen et al. [l( and Upadhyay et al. in 1998. Their 
works showed that determining the spin polarization at 
Fermi surface is essentially not an easy task. That leads 
to some definitions of spin polarization such as "tunnel- 
ing polarization" proposed by Tedrow and Meservey Q 
and "point-contact polarization" proposed by Soulen et 
al. [l(. One year later, Zhu et al. [J, Q and Kashi- 
waya et al. |6| have utilized the ideas to study the 
spin-polarized quasiparticle transport in ferromagnet/d- 
wave superconductor junctions. Zhu et al. 0, Q pre- 
dicted that conductance resonances occur in a normal- 
metal- ferromagnet/d- wave superconductor junction and 
in a following paper, they further studied the junctions 
by solving the Bogoliubov-de Gennes (BdG) equations 
within an extended Hubbard model which included the 
proximity effect, the spin-flip interfacial scattering at the 
interface, and the local magnetic moment. They have 
reported that the proximity can induce order parame- 
ter oscillation in the ferromagnetic region. In contrast, 
Kashiwaya et al. @ focused on the spin current and spin 
filtering effects at the magnetic interface. In the works of 
Zutic and Vails @, [1] , they first considered the effect of 
Fermi-wave vector mismatch (FWM) and have pointed 
out that if one neglects FWM, the effect of spin polar- 
ization invariably leads to the suppression of Andreev re- 
flection (AR). Among many other junction studies, Dong 
et al. Q studied a little different junction which forms a 
four layer sandwich, i.e., FM/I/d + is/d-wavc junctions, 
by taking into account the roughness of the interfacial 
barrier and broken time-reversal symmetry states. 



The pioneering works of Soulen et al. and Upad- 
hyay et al. have inspired several experimental studies 
El HI HI, H 0, EE M, 03 as welL Especially nor- 
mal and ferromagnetic metal/conventional superconduc- 
tor or s-wave superconductor (FM/s-wave SC) junctions 
have been intensely studied experimentally and theoret- 
ical modelings (Blonder-Tinkham-Klapwijk (BTK) for- 
mula [3 or its extension) had a good fitting with the 
conductance data. Recently Linder and Sudb0 [n| pre- 
sented a theoretical study of FM/s-wave SC junction 
that investigated the possibility of induced triplet pair- 
ing state in the ferromagnetic metal side. They have also 
used the BTK approach but allowed for arbitrary mag- 
netization strength and direction in the ferromagnet, ar- 
bitrary spin-active barrier, arbitrary FWM, and different 
effective masses in the two sides of the junction. As is 
expected, there is no retroreflection process when an ex- 
change field is present. However, they pointed out that 
retroreflection can occur under some conditions 

If one replaces the conventional superconductor by 
the high-temperature or <i-wave superconductor into the 
junction, it will occur several novel phenomena due to its 
d-wave pairing symmetry, complex band structure, and 
rich magnetic properties. Of particular interest, in the 
electron-doped side of cuprate superconductors (EDSC), 
it is strongly suggested that antiferromagnetic (AF) order 
may coexist with the d-wave superconducting order, es- 
pecially in the undcrdoped and optimally-doped regimes 
[20|. In this paper, we shall explore the possible novel 
phenomena in the FM/EDSC junction case, taking into 
account the interplay between antiferromagnetic order 
and spin polarization. The ideas and models developed 
in FM/CS junctions in the literature will be applied to 
the current FM/EDSC junction cases. 

This paper is organized as follows. In Sec. II, the 
basic formulation is given. We set up the condition of 
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the junction and generalize the BdG equations to in- 
clude AF order parameter. As the formal process, we 
utilize WKBJ approximation to obtain the more simple 
Andreev-like equations, which are then solved to deter- 
mine the four spin-dependent reflection coefficients (de- 
tailed derivations are given in Appendix [A"|) . Formulas 
of charge and spin conductances are derived . Sec. Ill 
are our main results and discussions. In Sec. III. A, the 
condition of midgap surface states was derived (details 
are given in Appendix [B]) . In Sec. III.B, the effect of 
FWM was studied. In Sees. III.C and III.D, we discuss 
the effects of spin-polarization and generalized effective 
barrier, respectively. It is shown that anomalous con- 
ductance enhancement can occur at low energies which 
could provide a further test for the existence of AF or- 
der in EDSC. In Sec. III.E, a more general formula for 
determining the spin polarization is proposed in terms of 
the experimental zero-bias conductance data. Finally in 
Sec. IV, a brief conclusion is given. 

II. FORMALISM 

Our formulation is given based on the following 
assumptions. We consider a point contact or pla- 
nar FM/I/EDSC junction where the superconduc- 
tor overlayer is coated with a clean, size-quantized, 
ferromagnetic-metal overlayer of thickness d, that is 
much shorter than the mean free path I of normal elec- 
trons. The interface is assumed to be perfectly flat and 
infinitely large. Considering I — > oo limit, the disconti- 
nuity of all parameters at the interface can be neglected, 
except for the SC order parameter to which the proximity 
effect is ignored [2l[. When SC and AF orders coexist, 
quasiparticle (QP) excitations of an inhomogeneous su- 
perconductor can have a coupled electron-hole character 
associated with the coupled k and k+Q [Q = (ir, n)] sub- 
spaces. QP states are thus governed by the generalized 
BdG equations [H H2 

Eui a (x) = H a u la (x) + J dx'A(s,r)v ls (x') + ^u 2lJ (x) 

Ev w {x) = J dxA*(s,r)u la (x)~H a v ls (x) + <S>v 2<f (x) 

Eu 2a (x) = $u la (x) + H a u 2a (x) - J dx'A(s,r)v 2s (x') 

Ev 2s {x) = §v 13 {x)- J dx'A*(s,r)u 2a (x') - H a v 2s {x), 

(1) 

where s = x - x', r = (x + x')/2, E\ o = -h 2 V 2 /2m - 
Ep' S —a J with J the exchange energy and a =\ (J.) for up 
(down) spin [a — —cr), and $ is the AF order parameter. 
A(s, r) is the Cooper pair order parameter in terms of rel- 
ative and center-of-mass coordinates. In the FM region, 
we define E% = h 2 q 2 F /2m = (h 2 q 2 F] /2m + h 2 q 2 F j2m)/2 
as the spin averaged value. It differs from the value in 
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FIG. 1: Schematic plot showing all possible reflection and 
transmission processes for an up-spin electron incident into 
the FM/I/EDSC junction. An AF order is assumed to 
exist in the EDSC. For convenience for a d-wave super- 
conductor, k x axis is chosen to be along the [110] direc- 
tion. The right-bottom inset shows a given Fermi wave 
vector = (k F ,k y ,k z ) and its coupled AF wave vector 
kf + Q = kfQ = (—k F ,k y ,k z ). Both vectors are tied to 
the Fermi surface, which is approximated by a square (thick 
line) . NR, AR, AF-NR, and AF-AR stand for normal reflec- 
tion, Andreev reflection, antiferromagnetic-normal reflection, 
and antiferromagnetic-Andreev reflection respectively. Their 
corresponding reflection angles are also shown. For the case 
of an incident down-spin electron, all spin indices just reverse. 

the superconductor, E F = h 2 k F /2m, to which a FWM 
can occur between the FM and EDSC regions Q. In JT]), 
the wave functions u\ and v\ are considered related to 
the k subspace, while u 2 and v 2 are related to the k + Q 
subspace. Comparing with the first and second lines of 
Eq. {T]), minus signs associated with the A(s,r) term in 
the third and fourth lines occur due to the symmetry 
requirement, A(k + Q) = — A(k), for a d 2 .2_ y 2-wave su- 
perconductor in k space. At Fermi level, the d a .2_ 1( 2-wave 

SC gap is A(k^) = Ao sin 29 with Ao the gap magnitude 
and 9 the azimuthal angle relative to the x-axis. 

In a d-wave superconductor, it's useful to consider a 
junction to which the superconductor surface is allied 
along the [110] direction. A thin insulating layer exists 
between the ferromagnetic metal and the superconduc- 
tor (see Fig. [T]) to which the barrier potential is assumed 
to take a delta function, V(x) — HS (x). Considering 
that an up-spin electron is injected into the FM/I/EDSC 
junction from the ferromagnetic metal side, there are 
four possible reflections as follows: (a) Normal reflection 
(NR): reflected as electrons, (b) Andreev reflection (AR): 
reflected as holes, due to electron and hole coupling in 
the k subspace. (c) Antifcrromagnctic-Normal reflection 
(AF-NR): reflected as electrons, due to the coupling of k 
and k + Q subspaccs. (d) Antifcrromagnctic- Andreev re- 
flection (AF-AR): reflected as holes, due to electron and 
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hole coupling in the k + Q subspace (see Fig. [T]). 

In addition to the effect of AF order, AR is largely 
modified due to the exchange field of ferromagnetic metal 
when electron is not normally incident into the EDSC 
region. Owing to the momentum conserved parallel to 
the interface, SnelPs law @, HH requires that 



q Fa sin 9 Na = q Fs sin 9 as = k F sin 9 Scr , 



(2) 



where On a, 9 as, and 9s a are the angles of NR, AR, and 
transmission into the SC respectively (see Fig. [T]). Inci- 
dent angle On a is typically not equal to the AR angle 6 as 
except when J = or for normal incidence. Assuming 
that there is no FWM and q F [ < k F < q F ^, ranges of 
six normal and Andreev reflection angles are < 0n\ < 
sin~ 1 (k F /q F f) = c2 , < A ] < sin' 1 (q F i/q F f) = 6 cl , 
and < 0s^,8si < sm~ (q F ^/k F ), while Qai and 0n[ 
can be any angles. For AF reflections, the angles = 
7r — Qa<j and Q N F a = On<t respectively. It is noted that 
when 6*tv| is within the range C \ < On^ < C 2, x com- 
ponent of the wave vector, J q 2 F ^ — k 2 F sin 2 9s-\ , becomes 

purely imaginary for the AR process [1, [f|. Although 
spin down electron as a propagating wave is impossible 
for AR, it can still transmit into the superconductor side. 

As emphasized by Kashiwaya et al. [6| , one can define 
two types of conductance in a FM, namely the charge and 
spin conductances. As a matter of fact, the normalized 
angle and spin dependent tunneling charge conductance 
is given by 



C qa = 1 — \R 
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a-a \Rac 



I n AF \ 
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where a\ = 1 and a-\ = L^X±-\ with Aif = 

cosOnt/ cosOsi, X21 = cos 0,4]./ cos 0sf, and L a = 



(q F /k F )(l — aJ/Ep). Detailed derivations of all four 

reflection coefficients (Rno-, Ras, Rn^i an d Ras) are 
given in Appendix [X] Similarly, the normalized angle 
and spin dependent spin conductance is given by 



C sa = l-\R 



Nn 



\Ras\ 



\R 



AF\ 
Ncr 



R 



AF I 



■(4) 



Comparing with the results of charge conductance in j3]) , 
due to the spin imbalance induced by the exchange field, 
different sign of Ras terms occurs in the spin conduc- 
tances. Consequently normalized total charge (spin) con- 
ductance is given by 

G q (s){E) = G q{s ^(E) ± G q ( a -)±(E), (5) 

where + (— ) sign is for charge (spin) channel and 

1 t (i 

Gq( s )a{E) — N / d9Na COS Ncr C q ( s y(E, 9 N a) Pa -(6) 

The lower and upper integration limits of a and are 
restricted by Sncll's law (as discussed before) or exper- 
imental setup. In practice, integration over two sepa- 
rate ranges of incident angle, i.e., < 1 9 Ncr | < c i and 



0c l < \Qn<j\ < 0c2 should be carried and results are added 
up to the total conductance. In ([6]), the normal-state 
charge (spin) conductance 



G 
where 



tt/2 

q(s)= I d9NaCos9 N AC m PT±CNiPi}, (7) 

-tt/2 



c 



\X X L C 



Na\ONcrj 



\l + X 1 L a +2iZ\ 



(8) 



with Z = mH/h 2 k F the barrier. In both © and 0, 
we have introduced a factor P a = (Ep + <jJ)/2Ep which 
can be interpreted as the probability of spin-er incident 
electron as a function of the exchange energy 0, H, . 
When J = 0, P T = P x = 1/2. 

In addition to the conductances, the normalized total 
charge (spin) current can be given by 



Iq(s) - ^(3)1 ± !q(s)l, 



(9) 



where 



iq(s)cr 



with 



1 fOO r[3 
—pj— / dE / d9No COS 9 NaC q ( s ) a {E 1 9NG-)PfTqF a 
q(s) J-oo J a 

(10) 



dE / d9 NrT cos9 N cr[CN]P}qF\ ± C N ]P]9f]]- 

-co J-tt/2 

(11) 



Charge and spin currents and their conversion are impor- 
tant probes for spin-related phenomena such as those in 
spin Hall effect. 



III. RESULTS AND DISCUSSIONS 

Both charge and spin conductances are important 
probes for tunneling in spin-polarized junctions. In this 
paper, we will focus on the charge conductance however. 
Moreover, for simplicity, all the results presented are for 
normal incidence (0^^ = 0). 



A. Midgap Surface States 

Detailed derivations of the midgap surface states 
(MSS) in the current FM/EDSC junction are give in Ap- 
pendix [B] Basically it is an extension of Hu's [2l| and 
Liu and Wu's [20| works. The boundary condition that 
leads to the MSS is the wave function ipNa(x = —d) = 
for a free boundary at x = —d. Consequently one obtains 
the following condition for the MSS (see Appendix [B]) : 



-2iki a d 



E, 



2$, 



(12) 
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where E± = E ± e' a with e' a = y/{E + a J) 2 - A 2 - $ 2 
incident spin-tr electron is assumed to have wave vector 
fcicr along the x direction. 

In case of J = 0, the result is reduced to our previous 
case without spin polarization [2(|. In case of J = $ = 0, 
the result is reduced to Hu's case 211, i.e., 



E 



E-e'' 



(13) 



where e' = \JE 2 — A 2 . The most crucial result of the 
above is that there exists a zero-energy state which is 
responsible for the ZBCP widely observed in hole-doped 
d x 2_ y 2-w&ve cuprate superconductors 21|. When J = 
but $ ^ 0, zero-energy state no longer exists such that 
the energy of the existing state is always finite (E = $ 
in the limit of d = 0). This leads to the splitting of the 
ZBCP. When J is also finite, there will be further effect 
caused by spin polarization although the splitting peak 
remains at E — $ in the limit of d — 0. It is interesting 
to note that beyond the quasiclassical approximation, a 
more accurate calculation for the surface bound-state en- 
ergies in d x 2_ y 2-w&ve and other unconventional cuprate 
superconductors was reported by Walker et al. j24| . 



B. Effect of Fermi- Wave- Vector Mismatch 

Tunneling conductances are in general strongly modi- 
fiedby the effect of Fermi- wave- vector mismatch (FWM) 
[1, [l!|. In our case, due to the presence of the AF or- 
der, the conductance spectra are somewhat different from 
those obtained by Zutic and Vails @ and Linder et al. 
[l9l ]. Here we introduce a parameter 



qt_ 



(14) 



to account for the effect of FWM. Both Lq greater and 
smaller than one cases are considered. As shown in 
Figs. [2]|4] for the normalized charge conductance G q , the 
effect of FWM is typically strong when Lq < 1, while 
it has relatively minor effect when Lq > 1 (That is, G q 
changes little from the no FWM Lq = 1 case.). As first 
pointed out by Blonder and Tinkham [2{| , FWM can be 
interpreted as a type of barrier which could enhance the 
conductance near zero bias. 

In our previous paper [2fj| . it was shown that ZBCP 
of a d x 2_ y 2-vjme superconductor can be split by the AF 
order No spin-active barrier 0, [lj|, external mag- 
netic field, and spin polarization effects were considered 
in our previous case though. Previously Zutic and Vails 
Q had given a detailed analysis of the FWM effect on 
the conductance in ferromagnet/s-wave and d-wave su- 
perconductor junctions. Here we show how FWM influ- 
ences the conductance in the current case and point out 
the key physics. Fig. [5] plots G q for various Lq with bar- 
rier Z — 0, AF order $ = 0.5Ao, and spin polarization 
X = 0.5 [see Eq. (fT7| for the definition of X], One sees 
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FIG. 2: Effect of FWM on normalized charge conductance 
spectra G q for various wave-vector mismatch value Lo with 
fixed barrier Z 
tion X = 0.5. 
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FIG. 3: Effect of FWM on normalized charge conductance 
spectra G q for various wave-vector mismatch value Lrj with 
fixed barrier Z = 1, AF order $ — 0, and spin polarization 
X — 0. This can be considered as the case of hole-doped 
cuprate superconductors without AF order and in the absence 
of spin polarization. 



that the effect of FWM is most noticeable at large FWM 
(Lq = 0.2 case) to which a ZBCP is developed, while 
the spectra are humdrum when Lq > 1. Since no barrier 
(Z = 0) is considered, no effect of AF order and spin 
polarization is seen in terms of peak splitting. Note that 
normalized zero-bias conductance is not equal to 2 due to 
the presence of AF order and spin polarization. In order 
to compare with the case of hole-doped high-T c supercon- 
ductors (without AF order), Fig. [3] plots G q at different 
values of Lq with $ = X = and Z = 1. One sees 
that ZBCP is largely enhanced by the FWM effect (see 
Lq = 0.2 case). Thus FWM can significantly enhance the 
number of midgap surface states near zero-bias voltage. 
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FIG. 4: Effect of FWM on normalized charge conductance 
spectra G q for various wave-vector mismatch value Lo with 
fixed barrier Z = 1, AF order $ = 0.5Ao, and spin polariza- 
tion X = 0.5. FWM causes the reduction of conductance at 
zero bias, while enhances the splitting peak associated with 
the AF order. 



Aiming to electron-doped cuprate superconductors, 
Fig. |4] shows the effect of FWM on the splitting peak 
when the AF order is present ($ = 0.5An). Here barrier 
Z = 1 and spin polarization X = 0.5. In contrast to the 
case of $ = in Fig. [3J FWM actually reduces the num- 
ber of midgap surface states near zero bias. At the same 
time, it enhances the strength of the splitting peak asso- 
ciated with the AF order. Following the idea of Blonder 
and Tinkham (2f| such that FWM can be interpreted as 
a type of barrier, the enhancement of ZBCP in Fig. [3j 
and the reduction of zero-bias conductance in Fig. 2] is a 
natural outcome at large FWM. 

In principle, the effect of FWM should be included 
when a serious calculation is performed for spin-polarized 
conductances. 



C. Effect of Spin Polarization 
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FIG. 5: Effect of spin polarization on normalized charge con- 
ductance spectra G q for various spin polarization value X 
with fixed barrier Z = 0, AF order $ = 0.5 An, and without 
FWM (L = 1). 
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FIG. 6: Effect of spin polarization on normalized charge con- 
ductance spectra G q for various spin polarization value X 
with fixed barrier Z = 1, AF order $ = 0, and without FWM 
(Lo = 1). This is considered an example of the hole-doped 
cuprate superconductor without AF order and FWM. 



In the literature, there exists different definitions of 
spin polarization. One example is the "tunneling polar- 
ization" proposed by Tedrow and Meservey [|[. In point 
contact experiment, the more suitable definition is the 
so-called "contact polarization" [l[ 



Pr 



N ] {E F )v F] ~N i {E F )v Fi 
Nf(E F )v F i + N l (E F )v Fl ' 



(15) 



where v Fa and N a (E F ) are respectively the Fermi ve- 
locity and DOS at Fermi level for spin-er electron. Since 
Lj cx N a (E F )v F<J , Eq. (fT5|) is identical to 



Pr 



h - h 



(16) 



However, the most natural definition of spin polarization 
is 



X 



N^(E F )-N i (E F ) 
N^E^ + N^EfY 



(17) 



In ballistic point contact situation, the electron density 
of states in the presence of an exchange field can be writ- 
ten as N a (E F ) = q FrT A/4:TT, where A is the area of the 
interface. Thus X = J I El with El = h 2 q 2 F /2m = 



{h 2 q 2 F] /2m + h 2 q F j2m)/2 ]l7j. In Sec. III.E, we will 
show that spin polarization X can be determined by a 
general formula in combination with the experimental 
conductance data. 

Note that current quasiparticle wave function of BdG 
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FIG. 7: Effect of spin polarization on normalized charge con- 
ductance spectra G q for various spin polarization value X 
with fixed barrier Z = 1, AF order $ = 0.5 An, and with- 
out FWM (Lo = 1). Low-energy anomalous conductance 
enhancement arises due to AF contributions (see text for de- 
tails). 



equations has four components which involve two com- 
ponents associated with the AF order. In the limit of 
Z = and without spin polarization (X = 0), nor- 
malized charge conductance has value 2 as expected (see 
Fig. [5]). With a finite AF order ($ = 0.5A ), the result- 
ing effective gap magnitude is about A ~ 1.12A (sec 
Fig. [5]). In general, at E < A, effect of spin polarization 
is to suppress the conductance. When FWM is absent 
(Lo = 1) together with Z = 0, normal reflection has no 
contribution and Andrcev reflection actually dominates 
the tunneling process for E < A [l8| . In our current case, 
Andrecv reflection involves contributions from both Ra 
and R^ F channels. 

The most interesting results occur when the barrier 
Z is finite. When the AF order $ = (as for the 
case of hole-doped cuprate superconductors) to which 
R^ F = Ra F = 0, ZBCP appears whose (normalized) 
strength is largely suppressed due to the strong spin po- 
larization effect (see Fig. [S]). However, as seen in Fig. 
when AF order is finite ($ = 0.5Arj), in contrast, the 
strengths of both the zero-bias conductance and the split- 
ting peak turn out to get enhanced by the strong spin 
polarization effect. This "anomalous conductance en- 
hancement" phenomenon is in drastic contrast as com- 
pared to the ZBCP associated with $ = case (Fig. [6]). 
These somewhat surprising results arise due to a signif- 
icant increase of |i?jy F | and at the same time, a signifi- 
cant decrease of \Ra\ I0r large X cases - a consequence 
of the interplay between AF order and spin polariza- 
tion. Since |i?)^ F | contributes positively to the conduc- 
tance, while |-R^ F | contributes negatively to the conduc- 
tance [see Eq. rcsultantly they cause the anomalous 
conductance enhancement at low energies (E < $). It 
should be emphasized that this low-energy conductance 



enhancement is not due to the spin-flip effect which is not 
considered in this paper. At higher energies, E > <!>, the 
conductances behave more normally such that they get 
suppressed due to the spin polarization effect. Anoma- 
lous conductance enhancement at low energies can serve 
as a test to see whether there is an significant AF order 
in electron-doped cuprate superconductors. 

Interface barrier and band structure are in general hav- 
ing strong effect on spin polarization. Kant et al. have 
built an "extended interface" model to illustrate the decay 
of spin polarization [12]. Besides, Mazin had a detailed 
discussion on the definition of spin polarization and band 
structure effects in spin polarization [26[. 



D. Effect of Effective Barrier 

In the study of the tunneling transition in Cu-Nb point 
contacts, Blonder and Tinkham [25( pointed out that bar- 
rier is not the only source for normal reflection and in a 
more realistic system, one should consider "impedance" 
or FWM as well which results in normal reflection even 
with no barrier present. They proposed an effective bar- 
rier Z e ff = [Z 2 + (1 — r) 2 /4r] 1 / 2 where r is the Fermi 
velocity ratio. They showed that effective barrier has 
an obvious effect on the conductance when E < Aq, as 
shown in Fig. 2 of Rcf. 



25] 



Here we generalize their idea 
to consider a spin, FWM, and angle dependent effective 
barrier Z e g @, [23 : 



[Z 2 + (l-L CT )74£a] 1/2 /cos% 



(18) 



where L a = qFal^F corresponds to the spin-dependent 
FWM. It is noted that we are not considering the spin- 
active barrier which has spin filtering effects and can 
lead to the ZBCP splitting [6|, [I2] . Instead we propose a 
possible alternative mechanism to account for the decay 
of spin polarization. Based on the generalized effective 
barrier, spin-up and -down particles experience different 
strength of effective barrier that causes spin-up and - 
down currents to decrease at different speed as compared 
to the current in the absence of barrier. Consequently, 
Z e g can modify the values of (and thus P c ) dramat- 

ically. With this strong effect at work, the decay of spin 
polarization should not be dominant by the spin-flitting 
process in the point contact spin polarization case. 

As seen in Eq. (fT5j) . Z e g can differ significantly from 
Z, especially when Z is small. Essentially their differ- 
ence can be measured by spin-polarized tunneling ex- 
periments. In Fig. [5J we compare the effects of Z and 
Z e g on the conductance with bare barrier Z set to zero 
and vary the FWM L value. For Z = and 6 Sa = 0, 
Z c fi = [(1- L a ) 2 /4L a ] 1/2 [see (fl8))]. In our case, we have 
also included AF order and spin polarization. The dif- 
ference is most noticeable when FWM is large (Lo = 0.2 
case). Since Z = 0, AF order and spin polarization have 
little effect at small FWM. However, when FWM is large, 
AF order and spin polarization can have a strong effect 
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0.5 1.0 1.5 2.0 

Normalized Energy E/A o 

FIG. 8: Effect of effective barrier Z e fj on normalized charge 
conductance spectra G q for various values of FWM Lq. AF 
order <£> = 0.5Ao and spin polarization X = 0.5. The bare 
barrier Z is set to zero, while Z a g is given by Eq. I|18|l. 



such that a splitting peak can develop at E m $ = 0.5Ao 
with the effective barrier Z e g (see Fig. [5]). This supports 
Blonder and Tinkham's idea of "impedance" mismatch 
which enhances the normal reflection. 



E. A General Formula for Determining the Spin 
Polarization 



Based on the phenomenon of Andrcev reflection. 
Soulen et al. [l[ proposed a formula for determining the 
point contact spin polarization P c [see Eqs. ([T5|) and (HU)] 
when the normalized zero-bias conductance data is com- 
pared. Their original form was 



G(0)/G N = 2(1 - P c ), 



(19) 



which is valid only when FWM is absent [llj. Since An- 
dreev reflection could be strongly modified due to the 
FWM effect, it's useful to replace Eq. (JEU) by 



G(0)/G N = 1 + \R A \ 



I R AF \ 



0--Pc), (20) 



where Ra and i?4 F are the AR and AF-AR coefficients 



respectively. Eq. (|20|) can be reduced to Eq. ([19)) when 
the exchange energy J is set to zero in Ra and the AF or- 
der $ is set to zero in Ra F - Note also that the parameter 
X should be set to zero when the "contact polarization" 
P c is determined under the idea of Soulen et al. 

Here we propose a more general formula for determin- 
ing the spin polarization: 



G{0)/G N = A< [ + A l , 



(21) 



where 



and 



Ai 



d0 Na cos Onv (1 + |.Rat I" 



i<T)n- 



(23) 



Here Ra<j 

TfAF 
H Aa 



R A a(L ,X,<t<,6 N(7 ) and R 



AF 
An 



(L ,X,$,6 Na ) with E = 0. Eq. (J2TJI is a natural 
result of our earlier formalism. It is regarded as the gen- 
eralization of Eq. (fTT))) of Soulen et al., which includes the 
effects of FWM, spin polarization, AF order, as well as 
the incident angle. 



IV. CONCLUSIONS 

Tunneling experiment provides a useful tool for prob- 
ing the properties of a superconductor such as the mag- 
nitude and symmetry of the superconducting order pa- 
rameter, quasiparticle density of states, and any existing 
competing orders. In fact, tunneling experiment is also 
a powerful probe for investigating the spin-charge sepa- 
ration in connection with the spin-injection techniques. 
This involves both charge imbalance and spin imbalance 
studies. 

In this paper, we have presented a detailed study of 
the tunneling conductance spectra of a ferromagnetic 
metal/electron-doped superconductor junctions, taking 
into account an AF order existing in the the electron- 
doped superconductor. Interesting result, such as low- 
energy anomalous conductance enhancement, occurs as 
a result of the interplay between AF order and spin po- 
larization (see Fig. [7]). These results in turn provide a 
further opportunity to test whether there is an significant 
AF order in electron-doped cuprate superconductors. 
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APPENDIX A: REFLECTION COEFFICIENTS 

Under the WKBJ approximation [U, EL EL HI, HI, 
HOj EL H2' the wave functions in the generalized BdG 
equations ([1]) can be approximated by 




\ 



(Al) 



A T = / ^0080^(1 + 011^^-^1^1^(22) 



Thus one obtains a set of Andrcev equations in the x 
direction, 



s 



H a u la (x) + A(k F )v la (x) + <&u 2a (x) 

A* (k F )u llT (x) - H a Vl 9 (x) + §V 2a {x) 
$>Ula(x) - H a U 2lJ {x) + A(k_FQ)u 2 tj (x) 
<&Vla(x) + A*(k FQ )u2cr(x) + H a V 2a (x), 

(A2) 



Ev ls {x) 
Eu 2cr (x) 
Ev 2s {x) 

where H a = — ltl kF 4- — a J and x is the coordinate nor- 

a m ax 

mal to the interface. The d x 2 _ y 2-wave SC gap A(kp) = 

— A(kpQ) = Ao sin 29 with 9 the azimuthal angle relative 
to the a;-axis. In obtaining Eq. (|A2[) , the Fourier trans- 
form of the Cooper pair order parameter A(s, r) from 
relative coordinate s to k space is assumed to take the 
form, A(k, r) = A(kp)8(x), with 8 (x) the Heaviside 
step function [2l|, [27} ■ 

Solving Eq. l[A"2]) . one obtains four eigenvectors which 
build up the spin-er wave function in the superconductor 
region (x > 0) @], 























+ c 2a 


;?)] 











C3c 




-ik x 



(A3) 



Here E± = E ± e a with e a = V ' E 2 - A 2 - $ 2 , A = 
A(kp), fc + = fc~ = kp cos9sa, and are coefficients of 
the corresponding waves. As pointed out by Blonder et 
al. [H], there is no need to normalize the coefficients as 
it just complicates the calculation. If we set $ = J = 
and normalize the coefficientSjit will reduce to the case 
for a typical N/I/S junction [U H3] . 

Since we consider that there is an AF order in the 
EDSC side, an incident electron from the FM side will 
have four possible reflections [2(| • The spin-cr wave func- 
tion in the FM side :( x < 0) with incident angle On a can 
thus be written as p, El, H3] 

iqw<j cos 8 Na x 



^Na{x) 



( 



\ 



RAse 

N(7 



iqF5 COS 9A5X 
iqFa cos 6 NlJ X 



1QF<t cos vn&x 



\ 



AF -iq FS cos 6 A 



n As e 



(A4) 



where Rno, Ras, Rnoi an d ^as are amplitudes of NR, 
AR, AF-NR, and AF-AR respectively. Applying the fol- 
lowing boundary conditions: 

IpNa (x) la—Q- = IpSa (x) | x=0 + (A5) 



4>Sa (x) | X =(H 



di>Sa (x) 



dlpNa- (x) 



\x=0~ 



2mH 

-p-VSaWl^ = dx 

the four reflection amplitudes (coefficients) are solved to 
be 

E—(l — L a \i a + 2iZg)B 



Rn<j — 



{l+L a X la + 2iZ g )D 



Ras 



r AF 
AF 



A(l + L s \ 2a +2iZ g )A 
(l + L a \ la +2iZ )D 

1 - L a \ la - 2iZ g 

1 + L a \ la + 2iZg 

A(l j- LjAicr — 2iZg)B 
(l + L s \ 2a ~2iZ e )D 

E„{l-L s \ 2s -2iZ e )A 



R 



(H 



L a X 2a - 2iZ 8 )D 



(A6) 



where 



.4 

B 
D 



= 2AL a X la [l - L a L s \ la \ 2s + AZ 2 

+ 2iZg(L a X la + L s \ 2a )] 

= 2L a \ la \2La\2aE + e{l + L%\l s )} 

= A 2 [(l — L a L s \\ a \ 2s + AZg) 2 

+ AZ 2 {L a X 1<T + L a \ 2s ) 2 ] 

+ [2L a X la E + AeZ 2 + e(l + Ll\ 2 1<7 )} 

x [2L a \ 2a E + 4eZ 2 + e(l + L 2 \ 2 2a )]. 



Moreover 



(A7) 
bar- 



Zg = Zj cos 9s a with the 
rier Z = mH /h 2 kp, \\ a = cosONa/ cos6sa> 
\ 2 s = cos 6 as I cos 6sa, and L a = 
\J(IfI^f — &{<If /kp){J I 'Epn). It is interesting to 
note in (|A6|) that R N ^ and arc proportional to the 
AF order as is expected. 



APPENDIX B: MIDGAP SURFACE STATES 



Following Ref . [2l| , we first assume that 



Via 



Ul<7 

Via 



(Bl) 



where 7 CT is the attenuation constant for ^(q^v)! < 
f $ 2 . With |B1]), Eq- GU) becomes 



|A(ki 




(h 
A 



A 


<I> 






-h 







vw ] 





-h 




i > 1 


<I> 


-A 




\ ^2<T / 



(B2) 



for the superconducting ovcrlaycr {x > 0). Here h = 
e'^ — oJ with e^. = ih 2 m~ l ^f a qp cos 9 The wave- vector 
components parallel to the interface are conserved for all 
possible processes. 

Solving Eq. (|B2p , one obtains double degenerate eigen- 
values E = A 2 + <i> 2 + e' a 2 — a J, where + (— ) corre- 
sponds to the electron- (hole-) like QP excitation. Simi- 
lar to the wave function (|A3|) . superposition of the four 
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eigenstates makes up the formal wave function for the 
superconductor overlayer (x > 0) 





( A ^ 
















o 


+ c 2(J 


:?)] 




^ $ ) 








e~ 1 " x e ik 



e~~ l " x e~ ik x 



(B3) 



Here E± = E ± e' a with e'„ = ^{E + aJ) 2 - A 2 - $ 2 , 
Cj are coefficients of the corresponding waves, and k + = 
k~ = kp cos 9 sa- At the interface, the wave functions of 
FM and superconductor meet ideal continuity ipN a (x = 
0) = 4>Sa(% = 0). After some algebra, the formal wave 
function for the FM overlayer is obtained to be (x < 0) : 











( e lkl ' x E+ \ 








e ~ik la .x 













+ C 2< y 










e tk la XQ J 




^ o J 





C3o 



I e ~^ la x E _ \ 
_ e ik lcr x ^ 
g ik\ a x ^jj 





e ik lt7 x E+ 





(B4) 



where it is assumed that incident spin-er electron has the 
wave vector k\ a along the x direction. Considering the 
free boundary at x = —d, ipN a (x = —d) = 0, one thus 
obtains the condition for the surface bound states: 



e~ 2iklad E + + e 2lklad E- = 2$. 



(B5) 
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